Proposal of a method for detecting dull images 
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Abstract - Arguments are proposed which show how images hardly percep- 
tible as hided in a clouding backgound can be revealed utilizing special patterns 
drawn by means of a simple mathematical procedure. Possible applications 
might be found in medical radiology. 
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Digitally recorded images often require to be processed so as to remedy lacks 
of readability due either to the limited capabilities of the recording devices or to 
the intrinsic features of the images itself, for instance, a poor contrast. Indeed, 
objects showing fading edges inevitably yield poor-contrast images, but also 
objects caracterized by sharp edges, if buried in a heavy background, cannot 
allow good-contrast images. Cases of this kind are common in medical radi- 
ology. Organs in the human body, such as the blood vessels or the digestive 
tracts, have indeed well-defined edges, but often X-ray opaque matters must be 
there injected to allow their X-ray imagery. Other cases of this kind happen in 
archaeology. Ancient masonries buried in the soil sometimes can be detected if 
the dim contrast they originate in the aerial pictures of archaeological sites is 
made out. 

Focusing our attention on the just outlined problems of imagery of hided- 
in-background sharp-edge objects, we consider the possibility of detect these 
images utilizing special patterns drawn by means of a mathematical procedure 
purposely devised jl. 

Utilizing a scale of greys or false colors, any monochromatic image lying on 
x, y plane is defined by a positive function <fi(x,y). When variation of <f> for 
variations of x, y in the range of the significant details of the image is very small 
with respect to actual values of </>, the image is dull, that is, barely perceptible. 
In practical cases, without significant limiting of generality, function <fi can be 
represented on x, y plane by means of contours, that is, lines where <p is constant. 
Let A be any contour and O any point on this contour. We chose point O as 
the origin of x, y axes and also as the origin of two orthogonal £, rj axes, the £ 
axis being tangent to line A. Transformation between these axes is 

£ = cos a ■ x + sin a ■ y, 

i] = — sin a ■ x + cos a ■ y, 

a standing for the angle between x and £ axes. By taking into account that at 
the point O we have d(j)/dS, = and putting dfyjdr) — S* , we get 
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2 ) Data on conventional methods to improve contrast in digital images can be found in 
Wikipcdia Ref. 1. 



dy d£ dy dr] dy 

which leads to 



S (x,y) = 




where 5 = \S*\. Since function S depends only on variation of on x,y plane, 
it allows to reveal the details of a dull image. To apply equation {T]) derivatives 
dip/dx and d<f>/dy must be evaluated numerically. Assuming, for instance, that 
is known in equally-spaced points with constant interval h (pixels), derivative 
in the central point among seven ones is 



4 = (-0i + 90 2 - 450 3 + 4505 - 90 6 + 7 ) /60h - /i 6 (7) /14fj, 

both on x and y axes. In the remainder term, 0^ 7 ^ means the seventh derivative 
of taken at some point interior to the x or y interval. Formulas for points out 
of centre are also available (Milne 1949). 

In order to show how this procedure works, let us consider the simple example 
of a disk of radius R in which value of function barely exceeds that in disk 
environment where only background exists. Accordingly, it is convenient to split 
function in two contributions, that is, 

<t> (x, y) = 0d (x, y) + 0s (x, y) . 

Assuming D be an arbitrarily-small positive constant, the disk spoken of can 
be represented by 

<1>d (x, y) = r/ ; = v t, (2) 

y/x^T7-R) /T D 



1 + cxp 

where T_d is a length related to the steepness of the disk border (0. Indeed, 
by letting r = yj ' x 1 + y 2 and assuming Tjj <C fl, in the disk outside, that is, 
for (r — R) /Yd ^ 0, we have <pu — and ~ 0s . On the disk contour, that 
is, for r = R, we have = D/2 + 0#- Like, in the disk inside, that is, for 
(r — R) /To <C 0, we have ~ D + 4>b- So, if 4>b is considered constant, it 
follows from {TJ and © that 

s (r) = D ex P [(r - R) /Tg] 

T D {l+exp[(r-R)/T D }} 2 ' 



3 ) For simplicity sake, in equation (f2j) we utilize a Fermi's function. Alternatively, the more 
steep Erf function could be considered. 



It is easy to check that 



D fd8\ fd 2 8\ D 



which mean that function 8 (r) shows a maximum on the disk contour. We have 
moreover 

D exp[(r-R)/T D ] _ D exp [(R - r) /Y D ] 



T D {1 + exp [(r - R) /T D }} 2 T D {l + ex p [(R - r) /T D }} 2 

which means that function 8 (r) is symmetrical across its maximum value 8 (R). 
Utilizing equation 

D cxp[a/T D ] 1 D 



T D {l + cxp[cr/r D }} 2 2 4IV 
the HWHM a of function 6 (r) is found to be 

a = 1.763 -T D . 

Consequently, it follows from that function <j>o (r) increases from 0.146 • D 
to 0.853 • D, for r increasing from R — a to R + a. 

The found results show that the disk <pjj (r) corresponds to a ring S (r) with 
the same radius R of the disk, the "thickness" a of this ring being related to 
the steepness of the disk border. It is important to point out that a result 
of this kind is fair for whatever figure on x, y plane more complex than a disk, 
provided that it can be defined by a contour line A allowing in all its points for 
a tangent £ axis. So, function 8 is expected to represent in general a contour 
pattern, with a "contour-thickness" a related to the sharpeness of the figure 
border. This result is based on the assumption of 0b constant. But even small 
variation of 4>b in i, y plane affect the actual value of 5 so that readability of 
the contour pattern might be worsed by a fragmentary backgound. A simple 
istance is 

1 B 

Bq and B\ standing for constant quantities and Tb for a length large with 
respect to Tb representing the variation of the image background along x-axis. 
Owing to (JTJ) , we have 

d^B Bi 

so that the condition for the contour readability obviously is 6 (R) > 8b, that 
is, remembering ([3]), 

D Tp 

4Bi - r B • 
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If this condition is fulfilled, contour patterns 8 reveal the presence of otherwise 
hardly perceptible dull images (0. This allows conclude that transformation {T]) 
could provide data sufficient for characterizing the objects under consideration 
even when their actual images are not available. In our opinion, the use of 
contour patterns might have some interest especially in medical radiology. 
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4 ) To apply transformation JT} when image is affected by a serious noise it could be 
convenient utilize previously a proper noise-filtering procedure (see Ref. 1). 
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